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ABSTRACT 

Recent observational data cast serious doubt on the widely held 
view that the sun's corona is heated by traveling waves (acoustic or 
magnetohydrodynamic) . u is here proposed that the energy responsible 
for heating the corona is derived from the free energy of the coronal 
magnetic field derived from motion of the "feet" of magnetic field lines 
in the photosphere. Stochastic motion of the feet of magnetic field lines 
leads, on the average, to a linear increase of magnetic free energy with 
time. This rate of energy input is calculated for a simple model of a 
single thin flux tube. The model appears to agree well with observational 
data if the magnetic flux originates in small regions of high magnetic 
field strength as proposed by Tarbell, Title and Schoolman. On combining 
this energy input with estimates of energy loss by radiation and of energy 
redistribution by thermal conduction, we obtain scaling laws for density 
and temperature in terms of length and coronal magnetic field strength. 


I. INTRODUCTION 


Following the general acceptance of the concept of a very hot corona 
in 1945 (Billings, 1966), Biermann (1948) and Schwarzschild (1948) inde- 
pendently suggested that the corona is heated by acoustic waves. However, 
recent observational data appear to make this hypothesis untenable. 

Athay and White (1979), in their analysis of UV spectroscopic data 

obtained by means of the 0S0-8 Spacecraft, argue that the acoustic wave 
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flux in the chromosphere cannot exceed about 10 erg cm s whereas the 
radiation losses of the transition region and corona require an energy 
flux of 10 erg cm s . Bruner (1980), from an independent analysis 
of 0S0-8 UV spectroscopic data, argues that most of the wave motion in the 
transition region is in the form of standing waves (probably evanescent 
waves) rather than traveling waves, so that the net acoustic flux is at 
least three orders of magnitude lower than that needed to heat "the corona. 

The presence of a magnetic field complicates the discussion. It is 
possible that the change is m'inor, in that energy still propagetes as 
waves, but these are magneto-acoustic waves and/or Alfven waves (Osterbrock, 
1961; Stein and Leibacher, 1974). If the waves are magnetoacoustic, we 
still have a problem in accounting for the required energy flux. If the 
waves are Alfven waves, then there is a problem in accounting for their 
dissipation (Stein and Leibacher, 1974) although Uchida and Kaburaki (1974) 
have argued that large-amplitude Alfven waves may convert to magneto- 
acoustic waves at coronal heights and thereby be dissipated. 

On the other hand, the close association between magnetic field 
strength and coronal heating suggests that the magnetic field may play 
an active role, rather than a passive role, in the energy transport 
process. For instance, it is notable that the Rosner-Tucker-Vaiana model 
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of hot coronal loops (Rosner, Tucker and Vaiana, 1978} seems to fit not 
only large-scale loop structures and active region loops, but also coronal 
loops produced by solar flares. Although this model deals only with the 
relationship between energy transport by heat conduction and by radiation, 
it does raise the question as to whether coronal loops are heated by some 
mechanism akin to that responsible for solar flares, which is believed to 
be the release of free magnetic energy by dissipative plasma processes 
(Sturrock, 1980). It is quite possible that the explosive nature of flares 
is due to the sudden rearrangement of magnetic field by an MHD instability: 
heating (and possibly acceleration) which occur during a flare may be due 
to dissipative processes which may also occur at a slower rate in the 
"quiet" solar atmosphere. 

It appears that there are two requirements for the propagation of 
preflare energy from the photosphere into the corona: there must be a 

magnetic field, and the footpoints of the magnetic field must be moved in 
such a way that the configuration is raised in energy from the current- 
free "ground state" to the current-carrying "excited state" which, in the 
early stage of excitation, may with good approximation be taken to be a 
force-free magnetic-field configuration. These requirements can be met 
equally well in the "quiet" solar atmosphere: the photosphere is always 

permeated by magnetic field, and the photosphere is always in motion due 
to granulation, supergranulation and other motions. The aim of this 
article is to estimate the resulting flux of energy from the photosphere 
into the corona and to incorporate this assumption into the Rosner-Tucker- 
Vaiana model so as to obtain estimates of the density and temperature of 
the plasma in a coronal loop in terms of the basic parameters of the loop. 

This possibility has recently been addressed also by Golub et al . 
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(1980). They begin with an analysis of observational data which indicates 

that there is a relationship between the plasma pressure in a coronal loop 

and the magnetic field strength in the loop, which they assert to be 
1 5 

p oc B ‘ . They estimate the rate of energy supply into a coronal loop in 
terms of the mean transverse (toroidal) magnetic field strength and the 
mean transverse (torsional) velocity at the photosphere,, but without con- 
sideration of the cause-and-effect relationship of these two quantities 
or their stochastic nature. Analysis of the data contained in their 
Fig. 2(b) yields a pressure-magnetic-field-strength relationship of the 
form 

p « B b , b = 0.77 ± 0.23. (1.1) 

This is a better fit to the relationship obtained in this article (b = 6/7, 
see Section 5), than the value b = 1 2/7 found by Golub et al . (1980). 
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II. TWISTED FLUX TUBE 


We wish to consider the response of a coronal flux tube to photos- 
pheric motion. Since the free energy available for heating the plasma is 
due to currents, we consider in particular a twisting motion such as 
would result from a rotation of either or both ends of the flux tube at 
the photosphere. For simpl it'.ty, we consider a tube which is thin compared 
with its length. We also simplify the calculation by considering a tube 
for which the central field line is a straight line and the tube has 
cylindrical symmetry. As indicated in Section I, we assume that the 
plasma density is sufficiently low that the field is approximately in a 
force-free state. 

It will be seen that the model we are considering is very close to 
that of Gold and Hoyle (1960) except that we are allowing for slow vari- 
ation of the field with respect to the z coordinate of cylindrical coordi- 
nates (z,r,<}>). On the other hand, we will simplify the model by consider- 
ing only the lowest order significant terms in polynomial expressions of 
quantities in terms of r. 

With these restrictions, we find that the solution of the equations 


Bx(V*x~B) =0 

may be expressed, to lowest significant order in r, as 

B z (z,r) = B q (z) - ^B; 1 (z) + b 2 B Q (z)Jr 2 , 
B r (z,r) = - j Bg (z) r, 

B (J) (z,r) = b B Q (z) r. 
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( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


We see that each field line is rotated about the axis by an angle Ax 
given by 

A x = b L, (2.4) 


when L (cm) is the length of the tube. 

If the radius of the flux tube is R(z), the magnetic flux passing 
through the tube is given by 

* " *B 0 R 2 - |B;' R 4 - |b 2 B 0 R 4 . (2.5) 


and the magnetic pressure at the surface of the tube is given by 


P 3 
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( 2 . 6 ) 


Let us now consider a flux tube, which is initially untwisted so 
that b = 0, immersed in a plasma so that the magnetic pressure at the 
surface of flux tube is balanced by an equal gas pressure outside the tube. 
Now suppose that the tube is twisted so that b f 0, but we require that the 
flux $ is unchanged and that the magnetic pressure p at r = R is unchanged, 
since it is still balanced by the same external gas pressure. Then we 
must expect that B (z) and R(z) will both change. It turns out, however, 
that to lowest order (quadratic) in b, R(z) is unchanged. On the other 
hand, B (z) is found to increase quadratically with b according to 

2 B°/3b 2 )b 0 = b o r2 - (2.7) 

The "free energy" AW of the flux tube is 'the magnetic energy of the 
twisted tube less the magnetic energy of the untwisted tube. .To lowest 
order in b, this is found to be 

AW = — fdz B 2 R 4 , (2.8) 

16 J 0 
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which shows that the free energy is distributed uniformly, along the length 
of the twisted flux tube. By using equations (2.4) and (2.5), we may re- 
express (2.8) as 


AW = 


$ 2 (AX ) 2 
16tt 2 L ’ 


(2.9) 
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III. STOCHASTIC MOTION 


We now consider the twisting of a flux tube due to stochastic motions 
at the photosphere. The general surface motion may be divided into two 
components, according to the equation 

v = - vy + Vx(An) (3.1) 

where n is the unit vector normal to the surface. The first term on the 
right-hand side denotes the "curl -free" component and the second term 
denotes the "vortical" component. 

We need to consider the transfer of energy from the photosphere into 
the magnetic field due to a random (horizontal) velocity field in the 
photosphere. A general treatment will be published at a later date. At 
this time, we undertake only a simplified analysis. 

We assume that the velocity field is statistically homogenous and iso- 
topic, and that half of the energy in the velocity field is in curl -free 
motion and half in vortical motion. Only the latter component contributes 
to twisting of magnetic field lines. We therefore consider an elementary 
contribution to the twisting by considering the rotation of one end of the 
flux tube where it meets the photosphere. Another contribution to the 
stored energy will come from the other end of the flux tube. 

If the vorticity at the center of the flux tube where it meets the 
photosphere is w, 

v^ = cor • (3.2) 

and the contribution to the twisting due to rotation at one end only is 
given by 

t 

A x = / co(t' )dt 1 . (3.3) 

o 
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This leads to the rate of change of'energy given by 

aw - it /l2d 2 \ 
dt ' ' At ' ' 


We see from equation (3.3) that 

t+At t+At 


<( (Ax) 2 ) ■ Jit’ Jit" <■>(’’ )»(t")> . 


t t 


This becomes 


oo 

((Ax 2 )) = AtJdxR(x) 


if 


(3.4) 


(3.5) 


(3.6) 


R(x) = <^co (t )co(t +x))> (3.7) 

and At is large compared with the range of x over which R(x) is significant. 
Hence if the "correlation time" x c is defined by 

oo 

JdT R(x ) = R(0)x c (3.8) 

o 

and if At » x , equation (3.6) becomes 

^iM!) 2 R(0) Tc - 2 < 0 ) 2 > x c . (3.9) 

If, considering only one end of the flux tube, the radius of the flux 
tube is R* and the field strength is B* at the photosphere, then 

) ' 2 "M 2 (<« 2 > v (3.io) 

The combination \co > R* is the mean square vortical velocity at the 
circumference of the flux tube, which is twice the mean square vortical 
component of the velocity over the flux tube. Hence it is the same as 
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the mean square velocity ( v ) , including both vortical and curl-free 
components. 

If we now consider both ends of the flux tube and assume that the 
field strength is the same at both ends, we find that the rate of increase 
of energy in the flux tube is given by 


* = !!i 

dt 4 ttL 



(3.11) 


Now suppose that R c is the mean radius of the flux tube in the corona 
defined by 

V = ttR 2 L (3.12) 


where V is the volume of the flux tube. Suppose also that the mean field 
strength B c is defined by 

3> = ttr2b c . . (3.13) 

_ O 1 

Then if (erg cm s ) is the mean rate of energy input into the flux 
tube, defined by 

(3.14) 

we see that 

£j= KB C L' 2 (3J5) 

where 

K = jLb.<v 2 >t. (3.16) 

4tt * c 
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IV. SCALING LAWS FOR CORONAL LOOPS 


We now consider a coronal loop and attempt to obtain scaling laws 

for the coronal density and temperature by investigating the energy 

-3 

balance. We dt by n Q (cm ) and T Q (K) the density and temperature at 
the top of the loop. For simplicity, we consider only loops which are 
sufficiently small that they are approximately isobaric. Then, if the 
mean density and temperature are n and T , 

V w 


We follow Rosner, Tucker and Vaiana (1978) in approximating the radi- 
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ation energy- loss function e R (erg cm s ) by 


e R ■ r "c T c' /2 


(4.2) 


wherein r ~ 10" 1 ®’®. 

If energy is deposited primarily near the top of a loop, it would be 

carried to the lower regions of the loop by heat conduction. If we denote 
-3 -1 

by (erg cm s ) the rate at which energy is extracted from the upper 
regions of the loop by heat conduction for transfer to the lower regions 
of the loop, then 


e H - 0 


T 7 / 2 L- 2 . 


(4.3) 


Noting that the "temperature scale height" is of order 1/2 L, and that 
the thermal conduction coefficient for a fully ionized plasma (Spitzer, 
1962) is approximately 10" 6 T 5 ^ 2 , we see that 0 * 10"^' 4 . 

If we now assume that 



(4.4) 
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we may obtain a relationship between n rt and T by assuming that e D « e u . 

oo K ri 

This relation is found to be 

T 0 - 0-V4 r l/4 s -5/8 „l/2 l 1/2 (4 . 5) 

Since the mean rate of energy input must be balanced by the mean 
rate of energy output (by radiation), we obtain what should be a more 
reliable relationship by equating the expressions in equations (3.15) 
and (4.2). This leads to the relation 

n l T ~c /2 = r_1 KB c L ~ 2 

for n c and T c - 

Alternatively, pursuing the assumption that « e^, we may set 
« £j and so obtain, from equation (3.15) and (4.3), the expression 

T . 0 -2/7 k 2/7 b 2/7 
0 c 


(4.7) 


V. DISCUSSION 


Equation (4.5)' has the same form as a relation obtained by Rosner, 
Tucker and Vaiana (1978). The more detailed model analysis of Vesecky, 
Antiochos and Underwood (1979) leads to a similar relationship with the 
numerical form 


V 
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•3.25 


nj / 2 L 1 ' 2 . 


(5.1) 


We find that this agrees with equation (4.5) if we adopt the proposed 

-0 15 

values of 0 and r and adopt 5 * 10 ’ i.e. T C /T Q =» 0.7, which seems 

not unreasonable. 

Equation (4.7) is of special interest in that it leads to a relation- 
ship between T Q and B c , independent of the length L. Although we have 
many magnetograph maps of the solar disk, unfortunately we do not know 
the magnetic field strength at the top of a coronal loop. Golub et al . 
(1979) have recently presented data for the coronal temperatures of thr«s_ 
typical features: x-ray bright points, active regions, and large-scale 
structures. The data presented by Golub et al . (1979) suggests the 
approximate numerical relationship 


T 


o 
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5.9 ± 0.1 



(5.2) 


In order for equation (4.7) to have this numerical form, we require that 
K * 10 15 - 2 ± We see from equation (3.16) that this requires the 

following combination of parameters at the photosphere: 


B*<v z > 


10 


16.3 ± 0.4 


(5.3) 


Tarbell, Title and Schoolman (1979) propose that most of the magnetic 
field in the solar atmosphere arises from small knots of intense field of 


* 


% 
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strength B* 3 1200 and that the r.m.s. velocity field in the photosphere 

* . 5 "1 

is about 10 cm s . Hence our theory leads to the empirical relation- 

i ship (5.2) if the correlation time for the photospheric velocity field 

has the value t c « 10 3,3 4 0,4 s, i.e. in the range 10 - 80 m. The 
lower limit of this range is comparable with current estimates (8 m) of 
the mean 'lifetime of granules (Allen, 1973). 

At this time, observational data concerning photospheric motions 
and the photospheric magnetic field, and observational data concerning 
magnetic-field strength at coronal levels, are insufficiently precise to 
provide a definitive check of the expression (4.7) for the coronal temper- 
ature in terms of the coronal magnetic fisld strength. In any case, 
relations derived in Sections IV and V depend on highly idealized assump- 

* tions concerning the structure and energy balance of coronal loops. The 
present estimates have been made simply to show that the heating mechanism 

t 

proposed in this article is not obviously ruled out by observational data. 

In order to obtain more detailed consequences for comparison with observa- 
tional data, we intend to incorporate the heating mechanism discussed in 
Section III into the model analyzed by Vesecky, Underwood and Antiochos 
(1979). In this way we shall avoid certain restrictive assumptions, such 
as constant pressure and the assumption that and we shall be able 

to include the important effect of magnetic-field geometry. 

Nevertheless, the fact that the models of Vesecky, Antiochos and 
Underwood (1979) lead to a scaling law similar to that derived, on a simpler 
basis, by Rosner, Tucker and Vaiana (1978) give some cause for optimism 

i 

that the scaling laws derived on the basis of the restrictive assumptions 

* of the current article may also prove to have wider applicability than 
might initially be expected. Ignoring the numerical coefficients, we 
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see from equations (4.6) and (4.7) that the density and temperature in 
a coronal loop will scale with the length of the loop and the coronal 
magnetic-field strength as 


so that 


b 4/7 l- 1 

c 

(5.4) 

r 2/7 
B c ’ 

(5.5) 

6/7 , -1 
c 

(5.6) 


It is notable that the index 6/7 (0.86) appearing in the scaling relation 
(5.6) is quite consistent with the value 0.77 ± 0.23 derived from the data 
of Golub et al. (1980), as noted in equation (1.1), and is in fact a 
better fit than the value 12/7 which they derived from their theoretical 
analysis. 

It will be interesting to see whether new data, such as may be obtained 
from the Solar Maximum Mission, will substantiate the above scaling laws. 
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